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Abstract 

On the basis of detailed ab initio studies the influence of strain on the anisotropy of the trans- 
port distribution of the thermoelectrics Bi2Te3 and Sb2Te3 was investigated. Both tellurides were 
studied in their own, as well as in their co-partners lattice structure to gain insight to the elec- 
trical transport in epitaxial heterostructures composed of both materials. It is shown, that the 
anisotropy of the transport distribution overestimates the experimental findings for Bi2Te3, imply- 
ing anisotropic scattering effects. An increase of the in-plane lattice constant leads to an enhance- 
ment of the transport anisotropy for p-doping, whereas the opposite occurs for n-doping. The 
recent findings and special features of the transport distribution are discussed in detail in relation 
to the topology of the band structures. 
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I. INTRODUCTION 



Thermoelectric (TE) materials have huge potential for power generation, heat pump- 
ing, and refrigeration. However, their practical application is restricted because of the low 
performance of the TE devices compared to traditional fossil fuel power generators and 
compressor-based refrigerators^. A high performance TE material has to be a good elec- 
trical conductor as well as a poor thermal conductor and, at the same time, possess a 
large Seebeck coefficient 3 -. Quantitatively the efficiency of TE devices is expressed by the 
dimensionless figure of merit ZT, 

ZT = a 2 aT/(K L + K e ), (1) 

where a, a, T are the Seebeck coefficient, electrical conductivity and absolute tempera- 
ture, and phonon and electron contributions to the total thermal conductivity, 
respectively. 

Bismuth and antimony tellurides, Bi2Te3, Sb2Te3, and the alloys based on these ma- 
terials play a significant role for thermoelectric technology. Already early studies of the 
(Bi 2 .,Sb2_ a; )Te3 compounds in the late 1950s reported ZT ~ 1 at room temperature^, con- 
firmed by further experiments 5 . This value remains, even to date, the maximal one available 
at room temperature for bulk materials. Current progress in nanostructure fabrication, in 
particular, epitaxial growth of high-quality superlattices^, encourages the possibility of sig- 
nificant performance improvement of TE devices. Venkatasubramanian et al. reported 
ZT ~ 2.4 and ZT ~ 1.45 at 300K for p-type and n-type Bi 2 Te3/Sb2Te3 superlattices, 
respectively^. 

These experimental advances motivated extensive theoretical studies of the electronic 
structure of the bulk bismuth and antimony telluride aimed to understand the possible 
origin of the increased thermoelectic performance in the multilayered structures. While 
in the previous years only few ab initio band-structure calculations of the bulk bismuth 
telluride^ could be mentioned, in the last decade various aspects of the electronic structure 
of both pure and doped bulk Bi 2 Te 3 and Sb 2 Te 3 as well as their transport properties were 
discussed in Ref. 10N22|. Ab initio studies of the electronic structure and the transport 
properties of Bi 2 Te3/Sb 2 Te3 superlattices were also reported 23 . 

An explanation of directional anisotropy of the transport properties in Bi2Te3/Sb2Te3 
superlattices could play a crucial role for the understanding of their increased figure of 



merit. Venkatasubramanian et al. 7 found a strong dependence of the anisotropy of the carrier 
mobility on both the superlattice period and the relative thickness of the constituents. The 
enhancement of the electrical conductivity parallel to the epitaxial growth direction, i.e. the 
trigonal axis of the rhombohedral lattice of bismuth and antimony tellurides, together with 
the possibility to suppress the lattice thermal conductivity Kl along this direction could 
provide the desirable ZT enhancement. 

In this study we concentrate on the anisotropy of the transport properties in the bulk 
Bi 2 Te3 and Sb2Te3 as a first step on a way of understanding the corresponding properties 
of the Bi 2 Te3/Sb 2 Te3 superlattices. Since epitaxial growth always implies lattice distortions 
we included, as discussed below, the effect of the lattice relaxation in our study. 

II. CRYSTAL STRUCTURE 

Both bismuth and antimony telluride possess the rhombohedral crystal structure with 
five atoms, i.e. one formula unit, per unit cell belonging to the space group D^ d (R3m). Ex- 
perimental lattice parameters^ are a r BiTe = 10.473A, 0BiTe = 24.17°, and a r s h bTe = 10.447A, 
OsbTe = 23.55°, were 9 is the angle between the rhombohedral basis vectors of the length a rh . 
In order to emphasize the layered character of this structure it is convenient to rearrange it 
into the hexagonal unit cell built up by three formula units, as shown in Fig{TJ The hexag- 
onal cell contains 15 atoms grouped in the 3 'quintuple' layers, Tel-Bi(Sb)-Te2-Bi(Sb)-Tel, 
where Tel and Te2 are two different crystal kinds of tellurium atoms. The 'hexagonal' lat- 
tice parameters are a h ^ Te = 4.384A, c h B e £ e = 30.487A, and a/gj£. e = 4.264A, c||f, e = 30.458A, 
for Bi2Te3 and Sb2Te3, respectively. In bismuth telluride the nearest interatomic distances 
between the individual monolayers inside the 'quintuple' blocks are 3.07A for Tel-Bi and 
3.25A for Bi-Te2. Two adjacent 'quintuple' layers in Bi2Te3 are separated by somewhat 
longer distance, 3.63A for Tel-Tel. In the antimony telluride these distances are 2.98A for 
Tel-Sb, 3.17A for Sb-Te2 inside the 'quintuple' blocks, and 3.74A for Tel-Tel between the 
blocks. 

In the Bi 2 Te3/Sb2Te3 multilayers atoms change their bulk positions due to the mismatch 
of the lattice parameters. The description of the realistic crystal structure of the multilayers 
is out of the scope of this study. Nonetheless, we modeled Bi2Te3 with the experimental 
lattice parameters and interatomic distances of Sb2Te3, and vice versa. Since both materials 
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FIG. 1: The rhombohedral unit cell superimposed with the hexagonal one to emphasize the layered 
character of the material. 

have very similar lattice parameter c hex along z-axis this variation is essentially compression 
and extension of the lattice in (xy) plane for bismuth and antimony telluride, respectively. 
We assume that one could estimate the effect of the lattice relaxation on the electronic and 
transport properties in the Bi2Tes/Sb2Te3 heterostructures from these two limiting cases. 



III. CALCULATIONAL DETAILS 



Calculations of the electronic structures were performed by means of the screened 
Korringa-Kohn-Rostoker Green's function method^ in the atomic sphere approximation 
(ASA) within the local density approximation of the density functional theory in the pa- 
rameterization of Vosko et al. 26 . It is generally recognized that the effects of spin-orbit 
coupling are mandatory for the correct treatment of the band structure in these materials. 
Therefore we used a fully relativistic version of the method based on the Dirac equation 27 . 
The obtained self-consistent band structures were used for the calculations of the transport 
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distribution a a p within the Boltzmann formalism assuming a constant relaxation time r— , 
^ (E) = r (2^S J dS ^)T> V(k) = - V i^(k), (2) 

J eJ(k)=S 

where £ J '(k) is the j-th band energy at the k-point of the Brillouin zone (BZ), a and {3 denote 
cartesian coordinates. We assume the relaxation time to be isotropic. In this approximation 
the transport anisotropy ratio o xx jo zz is independent on the relaxation time, and we do not 
have to specify it. 

The k-space integration over the isoenergetic surfaces was performed with the tetrahedron 
method on the Blochl mesh-2, of 96x96x96 points in the whole BZ. In the energy intervals 
of the width about 0.1 eV in the vicinity of both band edges the integration was refined by 
means of an adaptive mesh up to 4x4x4 times as dense as the original one. As a convergence 
criterium we used the values of the transport anisotropy ratio calculated from the effective 
mass approximation at the band extremal. 

IV. ELECTRONIC STRUCTURES AND TRANSPORT PROPERTIES 

The electronic band structures for Bi2Te3 and Sb2Te3 for both experimental and strained 
lattices are shown in Figj2j The positions of the high-symmetry points in the BZ of the 
rhombohedral structure are denoted in Figj3j Our results at the experimental lattice pa- 
rameters agree well with the previous ab initio studies of Mishra et air, Larson et alM for 
Bi 2 Te3, and of Eremeev et alM both for Bi 2 Te3 and Sb2Te3. At the same time the cal- 
culations made with the full-potential linearized augmented plane wave method (FLAPW) 
result in slightly different band structures for the bismuth and antimony tellurides at both 
the experimental^"— 1 ^ and strainecr^i lattices. The key question of the band structure 
of Bi 2 Te3 and Sb 2 Te3 is the position of the valence band maximum (VBM) and conduction 



band minimum (CBM) in the BZ. The calculations of Ref.|lll- ll4jj20| result in a six- valley 
both VBM and CBM located in the symmetry plane (TZU) in agreement with experiments 
for both the bismuth^ and antimony^ tellurides. Unlike these results, in our case the CBM 
of Bi 2 Te3, at both lattice parameters, asiTe^ and asbTe, is a two-valley minimum located 
on the symmetry line TZ, similarly to Ref.|9j. For Sb 2 Te3 at the experimental lattice pa- 
rameters we found a direct band gap located at the center of the BZ, while at the larger 
in-plane lattice constant both six- valley VBM and CBM lie in the symmetry plane (rZU). 



Bi 2 Te 3 energy band structures 



Sb 2 Te 3 energy band structures 
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FIG. 2: Band structures of Bi2Te3 and Sb2Te3 along symmetry lines for both experimental (solid) 
and strained (dashed) lattices. Energies are given relative to the VBM. 
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FIG. 3: Brillouin zone of the rhombohedral lattice. 

In contrast to our results Thonhauser— found that the increase of the lattice parameters in 
Sb2Te3 led to the formation of a direct band gap at the T point. On the other hand, the 
negative hydrostatic pressure discussed in Ref. [15[ implies an increase of both in-plane and 
out-of-plane lattice parameters in comparable degree, while in our study essentially the first 



one is included. Additionally, calculations in Ref. [15j were performed for optimized atomic 
positions with respect to the total energy, which can affect the band structure of antimony 
telluride^. 

As already discussed-' 11 , these differences in the band structures are probably due to the 
non-spherical part of the potential, which is not included in the ASA. At the same time, as 
discussed below, these differences have no significant impact on the transport distribution 
a(E). The details of the band structures for all four systems are compiled in Table I. 

The calculated transport distribution of Bi 2 Te3 at both the experimental, a = ciBiTe, and 
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TABLE I: Band structure parameters: Band gap in eV, positions of VBM and CBM in crys- 
tallographic coordinates, effective masses in electron mass units, principal axes e, in cartesian 
coordinates, and transport anisotropy ratio from the effective mass approximation. 



Bi 2 Te 3 


Lattice parameters 


a = &BiTe 


a = 


O-SbTe 


Gap(eV) 


0.105 


0.129 


Extremum 


VBM 


CBM 


VBM 


CBM 


Position 


0.517 0.366 0.366 


0.173 0.173 0.173 


0.405 0.405 0.335 


0.151 0.151 0.151 


Effective masses 










Tfl\ 


-0.024 


0.178 


-0.039 


0.154 


m 2 


-0.134 


0.178 


-0.077 


0.154 


m 3 


-1.921 


0.835 


-0.207 


1.370 


Principal axes 










ei 


0.500 -0.867 0.000 


1.000 0.000 0.000 


0.866 0.499 -0.024 


1.000 0.000 0.000 


e 2 


0.600 0.346 0.721 


0.000 1.000 0.000 


0.500 -0.867 0.000 


0.000 1.000 0.000 


e3 


0.625 0.361 -0.693 


0.000 0.000 1.000 


0.021 0.012 0.999 


0.000 0.000 1.000 




5.452 


4.700 


4.020 


9.013 


Sb 2 Te 3 


Latt. param. 


a = asbTe 


a = 


a-BiTe 


Gap(eV) 


0.090 


0.140 


Extremum 


VBM 


CBM 


VBM 


CBM 


Position 


0.000 0.000 0.000 


0.000 0.000 0.000 


0.547 0.392 0.392 


0.004 0.020 0.020 


Effective masses 










mi 


-0.054 


0.045 


-0.039 


1.124 


m 2 


-0.054 


0.045 


-0.083 


1.774 


m 3 


-0.102 


0.114 


-2.046 


6.861 


Principal axes 










ei 


1.000 0.000 0.000 


1.000 0.000 0.000 


0.500 -0.867 0.000 


-0.316 -0.183 0.931 


e 2 


0.000 1.000 0.000 


0.000 1.000 0.000 


0.594 0.343 0.727 


0.500 -0.867 0.000 


e3 


0.000 0.000 1.000 


0.000 0.000 1.000 


0.630 0.363 -0.686 


0.806 0.465 0.365 


&xx/ ' zz 


1.889 


2.507 


2.397 


2.080 
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FIG. 4: Transport distribution cr xx (E) and a zz (E) and the transport anisotropy cr xx /cr zz for Bi2Te3 
at the experimental and compressed lattice parameters. Crosses at the band edges mark the o xx ja zz 
ratio derived from the effective mass model^ using the parameters of Table I. 

compressed, a = asbTe, lattice parameters are shown in FigH] together with the anisotropy 
ratio a xx /a zz . In terms of the rigid band model the energies below and above the band gap 
simulate p- and n-doping respectively. While for p-doping close the VBM a xx /a zz varies 
smoothly approaching the limiting value, the ratio increases drastically and forms a promi- 
nent two-peak structure for n-doping case. This structure originates from the two topological 
transformations of the constant energy surfaces in the conduction band. FigJS^a) shows the 
contour plot of e(k) for Bi2Te3 at a = asiTe in the plane (rZU) for energies to 0.19eV 
relative to the conduction band edge. The main features of the band structure are the global 
CBM on the line TZ, the local conduction band minimum (LCBM) at (0.666,0.602,0.602), 
and two saddle points si at (0.722 ,0.667, 0.667) and S2 at (0.493,0.461,0.461), in crystallo- 
graphic coordinates. The saddle point si occurs at E — E C bm = 0.04eV and causes the first 
peak of (T xx /a zz , while s 2 appears at 0.17eV and forms the second peak. At room temper- 
ature the chemical potential would be located in the saddle point s\ or s 2) for an electron 
carrier concentration of about N = 3.0 x 10 19 cm -3 or iV = 1.5 x 10 20 cm -3 , respectively. 
The location of the LCBM in our band structure calculation is close to the position of the 
six- valley CBM reported in Ref. 11 1^ . I9IJ20I ] . and the energy difference between these two 



minima, Elcbm — EcBM=0.04eV, is quite small. At the same time the saddle point S2 close 
to the local band maximum at Z was found as well in agreement with these calculations. A 
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FIG. 5: Contour plots of e(k) at the experimental lattice constants in the plane (TZU) 

(a) Bi2Te3, 10 isolines for (E — Ecbai) at to 0.19eV with a constant increment (dotted), addi- 
tionally, 2 isolines at E — Eqbm = 0.04eV and E — Eqbm = 0.17eV with the saddle points s\ and 
S2, respectively (bold), the positions of the CBM and the LCBM are marked with crosses. 

(b) Sb2Te3, 10 isolines for {E — Evbm) at -0.24 eV to with a constant increment (dotted), 
additionally, 1 isoline with the saddle point s (bold), the position of the VBM is marked with a 
cross. 

slight shift of the energies at the CBM and the LCBM would bring these two band structures 
in accordance. Since the saddle point s\ lies close to the line connecting the CBM and the 
LCBM this modification would not affect remarkably the band structure topology. More- 
over, the transport anisotropy cr xx /a zz at the LCBM from the effective mass approximation 
is 4.95, which is fairly close to 4.7 at the CBM. This indicates that the transport anisotropy 
of the bismuth telluride at the experimental lattice parameters is stable with respect to small 
modifications of the band structures with a two- and six-valley CBM, respectively. 

In Bi2Te3 the in-plane compression of the lattice parameters from a = asiTe to a = asbTe 
increases the transport anisotropy ratio at the conduction band edge remarkably. Within 
the effective mass approximation this can be explained by the enhancement of the ratio 
Tn±/m\\ due to the expansion of the BZ in xy plane. At the same time a xx /a zz decreases at 
the valence band edge. In this case the compression of the lattice results in a re-orientation 
of the longest axis of the effective mass ellipsoid near the VBM closer to the z axis, which 
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FIG. 6: Transport distribution and the transport anisotropy for Sb2Te3 at the experimental and 
expanded lattice parameters. Crosses at the band edges mark the o xx ja zz derived from the effective 
mass model^ using the parameters of Table I. 

lead to an increase of the transport anisotropy ratio, and at the same time decreases the 
anisotropy of the dominating effective masses m 2 and m 3 (see Table I). The last effect prevail 
over the enhancement of u xx /u zz due to the rotation of the effective mass ellipsoid. Figj6] 
shows the transport distribution and the anisotropy ratio of Sb2Te3 at both a = asbTe and 
Q> — O'BiTe- The kink of the anisotropy ratio at the experimental lattice parameters is induced 
by the saddle point s =(0.831, 0.784, 0.784) at E — Evbm = 0.116 eV, which corresponds to 
a hole carrier concentration of N = 5.8 x 10 19 cm~ 3 at 300K. This topology is illustrated in 
Figj5]^b), which shows the contour plot of e(k) in the valence band of the antimony telluride 
for energies -0.24 eV to relative to Evbm- The in-plane expansion of the lattice parameter 
increases the density of the occupied states near the valence band edge and suppresses the 
kink. The transport anisotropy ratio increases with the in-plane compression of the BZ at 
the valence band edge, and, at the same time, decreases at the conduction band edge due 
to the larger angle between the z axis and the largest axis of the effective mass ellipsoid, 
similarly to the discussed changes in bismuth telluride. 
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V. CONCLUSIONS 



On the basis of ab initio electronic structures obtained with the fully relativistic KKR 
method we studied the anisotropy of the transport properties of the bismuth and antimony 
tellurides in the constant relaxation time approximation within the Boltzmann formalism. 
In addition to the systems with the experimental lattice parameters we modeled bismuth 
telluride within the lattice of Sb2Te3, and vice versa. We found that a decrease of the 
in-plane lattice parameters increases the transport anisotropy for the n-doping and, at the 
same time, decreases the anisotropy for the p-doped case. This effect can be understood 
within the effective mass approximation at the valence band maximum and conduction band 
minimum, respectively. 
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